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Ellipsoidal micron-sized colloidal particles can oscillate spontaneously when trapped in a focused 
laser beam. If two oscillating particles are held in proximity their oscillations synchronize through 
hydrodynamic interactions. The degree of synchronization depends on the distance between the 
oscillators and on their orientation. Due to the anisotropic nature of hydrodynamic coupling the 
synchronization is strongest when particles are arranged along the direction of oscillations. Similar 
behavior is observed for many oscillating particles arranged in a row. Experimental observations 
are well reproduced with a model that uses a phenomenological description of the optical force and 
hydrodynamic interactions. Our results show that oscillating ellipsoidal particles can serve as a 
model system for studying hydrodynamic synchronization between biological cilia. 
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Synchronization appears at all length scales from 
atoms to macroscopic bodies and is ubiquitous in living 
systems [1]. For example at microscale it plays an impor¬ 
tant role in the motion of microswimmers such as proto¬ 
zoa [2], algae Q or spermatozoa Q which use flagella or 
cilia for their propulsion. A green alga Chlamydomonas 
swims with two flagella which it moves in a coordinated 
way reminiscent of a human doing the breaststroke. The 
coordination of flagella is crucial for swimming along a 
straight line [3]. Ciliates like Paramecium coordinate the 
beating of their cilia to form metachronal waves which 
greatly enhance their swimming efficiency [E[. Similar 
waves are also observed in cilia that cover respiratory 
epithelia and clear mucus from the airways [ij- Hydro- 
dynamic synchronization can even be observed between 
flagella of two separate sperm cells swimming beside each 
other Q. In large numbers hydrodynamically synchro¬ 
nized spermatozoa can form intriguing vortex patterns 
[ 3 ] • The formation of metachronal waves has been recre¬ 
ated in an artificial system consisting of microtubule bun¬ 
dles and kinesin motors [8]. 

Synchronization of nearby swimming microorganisms 
with waving tails was first studied by G. Taylor who 
showed that it can lead to a reduced dissipation @. How¬ 
ever, this finding does not explain the kinematic mecha¬ 
nism that keeps the tails synchronized. A difficulty lies in 
the temporal reversibility of the Stokesian hydrodynam¬ 
ics while synchronization is irreversible by nature na. 
This can be overcome either by breaking the temporal 
symmetry in the driving mechanism [lj], [l2| or by de¬ 
scribing a cilium with more than one degree of freedom 
pj . In addition, synchronization can also result from 
indirect coupling between cilia caused by a rocking mo¬ 
tion of the cell body Bill , or from inertial effects at 
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non-zero Reynolds number 0. 

Hydrodynamic synchronization was studied in model 
system containing spherical colloidal particles actively 
driven along closely spaced circular trajectories via feed¬ 
back controlled optical tweezers [17-21J. Cicuta and co¬ 
workers demonstrated an in-phase and anti-phase syn¬ 
chronization where the transition between the two modes 
was controlled by changing the shape of the driving po¬ 
tential 0. In all these experiments the driving mecha¬ 
nism required video particle tracking and computer con¬ 
trolled feedback. 

Alternatively hydrodynamic synchronization can be 
studied in a system of optically trapped non-spherical 
particles. A recent study reported synchronous rotation 
of two particles in vortex beams [22|. Elongated parti¬ 
cles in a focused beam can also exert oscillatory motion 
[23l HU. Their oscillations are induced by a combination 
of the gradient force that pulls the particle towards the 
beam center and a non-conservative force due to the radi¬ 
ation pressure. In this Letter we report on experiments in 
which such oscillating ellipsoidal particles were exploited 
to study the hydrodynamic synchronization between two 
or more closely spaced autonomous oscillators. 

Experiment. Ellipsoidal prism shaped particles 
(Fig.[]Ji) were fabricated from a layer of photoresist (SU- 
8, MicroChemicals) by maskless photo-lithography using 
direct laser structuring (LPKF ProtoLaserD) J25j. Par¬ 
ticles were dispersed in water and sealed in a thick sam¬ 
ple cell. The experiment was conducted on a setup built 
around an inverted microscope (Zeiss Axiovert) equipped 
with a multi-trap laser tweezers system (Aresis Tweez) 
using 1064 nm IR laser (Coherent Compass) and high nu¬ 
merical aperture water immersion objective (Zeiss Achro- 
plan 63 x 0.9 NA) [26]. Kinematic information on particle 
motion was obtained from video recordings captured with 
a CMOS camera (Pixelink PL-B 741) and subsequently 
analyzed by proprietary particle tracking software |27j . 
At the beginning of each measurement one or multiple 
particles sedimented at the bottom of the experimental 
cell were trapped and pushed against the upper cell wall. 
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FIG. 1: a) SEM image of a particle with major radius 8 p m, 
minor radius 2.5 pm and thickness 2.2 pm. b) p denotes the 
angle between the particle’s long axis and the z-axis, x the 
distance between its center of reaction and the center of the 
laser beam. 


In order to control the orientation of a particle each trap 
consisted of a principal beam and a secondary beam with 
20% of the main beam’s power that was shifted against 
the main beam by 2 p m perpendicular to the x-axis. 

Single particle oscillations. Once a trapped particle is 
pushed against the upper cell wall (Fig. [TJd) it reaches a 
steady state in which it persistently oscillates about the 
center of the trapping laser beam. The oscillations take 
place in a plane orthogonal to the flat surface of the parti¬ 
cle. Figured shows a time series of images demonstrat¬ 
ing the oscillatory motion. The position of the particle 
center as a function of time is shown in Fig. [2t>- Since all 
optical forces are proportional to the light intensity we 
expected and found a linear dependence of the oscillation 
frequency on the laser power (Fig. [2t). 

The dynamics of an ellipsoidal particle in the optical 
tweezers beam can be described as follows. If we define 
the particle position as its center of reaction [28] its mo¬ 
bility tensor becomes diagonal (a force causes no rotation 
and a torque no translation) and the equations of motion 
are 


x = MttF , p = M rr T (1) 

where p is the angle between the particle’s long axis and 
the z -axis and x the distance between the center of reac¬ 
tion and the axis of the laser beam (Fig.[fl)). We propose 
the following expressions for the force and the torque act¬ 
ing on the particle 

F - —Ax + Bp 

( 2 ) 

T = -Cx + Dp — Ep 3 . v ' 

The coefficient A describes the restoring potential of the 
laser trap, B is the effect of tilt - the radiation pres¬ 
sure on a tilted particle exerts a lateral force; C is the 
restoring torque of the trap boundary; D describes the 
torque on a particle pressed against a wall that causes an 
instability and E is a restoring torque that counteracts 
the instability. The equations of motion © can then be 



FIG. 2: a) Time sequence of an oscillatory motion of an el¬ 
lipsoidal particle trapped by laser beam, b) Dynamics of the 
oscillating particle in y direction, c) Particle’s oscillation fre¬ 
quency as a function of laser power. The solid line represents 
a linear fit to the data, d) Oscillatory period of an ellipsoidal 
particle comparing the experimental and theoretical data for 
a = 0.8 and /3 = 1.5. e) Stable limit cycle for same values of 
a and /3 as in d). 


written in the form 

x = 0 (—ax + ep) 

p — 0 (—x + /3p — Sp 3 ) . 
e 

Of the five parameters O determines the oscillation fre¬ 
quency, S the amplitude and e the amplitude ratio be¬ 
tween x and p. The dynamical behavior of the system 
is governed by the two essential parameters {a and /?). 
The dynamical system described by Eq. [3] always has a 
fixed point at (0, 0) and for a/3 > 1 it has two more at 
= hy/(af3 — 1 )/(a5)(e/a, 1). The fixed point at (0,0) is 
stable if a > (3 and a/3 < 1 . At a = (3 it undergoes 
a Hopf bifurcation and becomes unstable, encircled by a 
limit cycle. At a/3 = 1 there is a pitchfork bifurcation 
(transition from 1 to 3 fixed points) and for a/3 > 1 (0,0) 
becomes a saddle point. The remaining two fixed points 
are stable if <a+2/3—3/a > 0 and unstable otherwise. The 
latter transition is a subcritical Hopf bifurcation. The 
dynamical regimes of the system are shown in Fig. [3l 
Figure [2]T shows a measured oscillation cycle obtained 
by averaging the experimentally obtained particle posi¬ 
tion as a function of time over many periods. The shape 
of this oscillation cycle is visibly asymmetric and com¬ 
pares well to the trajectories obtained from model equa¬ 
tions in the region with 3 fixed points. Equations of mo- 
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FIG. 3: Bifurcation diagram for a single particle described 
by Eqs. [3] Depending on the parameters a and /?, the model 
can exhibit one or three fixed points and a limit cycle that 
describes autonomous oscillations. In addition, the diagram 
shows (dashed line) the regions in which two particles show 
in-phase and anti-phase synchronization for parameter values 
corresponding to longitudinal arrangement at large distances 
(■ crr/ctt = 0.5, ctr/(cctt) = 0.5, crtc/ctt = 1.5). 


tion © were solved numerically for a = 0.8 and /3 = 1.5. 
The limit cycle of the oscillator for the same parameters 
is shown in Fig. [2k- From this comparison we conclude 
that an oscillating particle can be adequately described 
by the phenomenological model we introduced here. Fur¬ 
thermore, it has been demonstrated that an oscillator 
with two degrees of freedom is a good candidate for hy¬ 
drodynamic synchronization in fl3| . 

Two particle synchronization. To study hydrodynami- 
cally mediated synchronization between independent os¬ 
cillators we trapped two ellipsoidal particles in separate 
beams and brought them in close proximity. Measure¬ 
ments were performed for three different alignments of 
the particles with respect to the direction of oscillation: 
longitudinal (inset Fig. [4 Jd) , diagonal (Fig. 0k) and lateral 
alignment (Fig.0)i). An example of particle positions vs. 
time is shown in Fig. 0ki for an interparticle distance of 
8 /am (half the particle size) and longitudinal configura¬ 
tion. Synchronized oscillation is clearly visible. 

The degree of synchronization can be visualized in 
a scatter plot of particle positions x\ vs. x 2 (Fig. 0 Jd- 
d). For an interparticle distance of 15 /am the graph 
shows a strongly correlated motion for longitudinal align¬ 
ment (Fig. 0k>) indicating that oscillating particles syn¬ 
chronized in-phase. For diagonal alignment the motion 
is still synchronous although the correlation is weaker 
(Fig. 0k). For the lateral alignment the motion is not 
correlated showing that the oscillations are out of syn¬ 
chrony (Fig. 0]1). This observation is in accordance with 
the anisotropy of Oseen’s tensor which is even more pro¬ 
nounced in the presence of a no-slip boundary [29). We 
quantitatively evaluate the degree of synchrony by cal¬ 
culating the Pearson correlation coefficient as a function 
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FIG. 4: a) Position oscillations of two particles in longitudinal 
arrangement (inset in b)) at a distance of 8 /im, laser power 
56 mW. b)-d) Synchronization strength at distance 15 /am for 
three arrangements: b) longitudinal, c) diagonal, d) lateral. 
The insets show the position of the particles. Laser power: 
56 mW. e)-f) Correlation coefficient r x as a function of the 
interparticle distance d for different arrangements showing 
(e) in-phase synchronization (laser power 28 mW) and (f) 
anti-phase synchronization (laser power 42 mW). g)-h) The¬ 
oretically calculated correlation coefficients with parameters 
a — 0.8, (3 — 1.5, e = 14 /im (g) and e = 6 fim (h). 


of interparticle distance (Fig. 0k). In one instance the 
particles also showed robust antiphase synchronization 
(Fig. 0f). The type of synchronization was largely in¬ 
variant of the laser power. We conclude that the tran¬ 
sition from in-phase to anti-phase synchronization was 
governed by the detailed properties of the laser beam in 
the vicinity of the surface. 

To understand the hydrodynamic synchronization we 
extend the phenomenological model © by taking into ac¬ 
count the fluid-mediated interaction. The force or torque 
which sets one particle and the surrounding fluid in mo¬ 
tion also has an effect on the velocity and angular velocity 
of the second particle. Their equations of motion can be 
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written as 

£1,2 = MttF 1,2 + Ctt^2,i + CtrT2 : i 
01,2 = MrrTi,2 + C#T ^ 2,1 + CrrT 2: 1 

where x$ and pi denote the deflection of i-th (i = 1,2) 
particle, F$ the force and Ti the torque acting on it. Ctt, 
Ctr, • • • denote the coupling coefficients (off-diagonal 
elements of the two-particle grand mobility matrix). By 
introducing the reduced coefficients ctt = Ctt/Mtt , 
CRR = Crr/Mrr , CTi ? = Ctr/Mrr , CRT = Crt/Mrt 
and inserting F and T from Eq. [2] we obtain 


£1 = Q x (-axi + e<^i) + c T tH 2 (-<^2 + e^ 2 ) 

+ ct#H 2 (—-x 2 + Pp2 — bp 2) 

01 = fii(—xi + flip I - dpi) + crt^2{-0'X 2 + e<p 2 ) 

+ Ci?i?^2(--^2 + PP2 - dpi) 

(5) 

for particle 1 and analogous equations for particle 2. In 
this expression we assumed that the two oscillators differ 
only in their characteristic frequencies (11 1 and H 2 ). A 
numerical solution of the equations of motion (|5]) for two 
identical oscillators (Hi = fi 2 ) reveals that if the particles 
are only coupled through the crt coefficient, they will 
synchronize in anti-phase except for a narrow region with 
small a. The coefficients crr and crt always lead to in- 
phase synchronization and crr always to anti-phase. As 
a consequence, the nature of synchronization will depend 
sensitively on the ratio between coupling coefficients. 

We determined the coupling coefficients numerically by 
modeling the particles as ellipsoids orthogonal to a wall. 
We solved the hydrodynamic problem with a boundary 
element method (BEM) by representing each particle 
with 512 elements and incorporating the no-slip bound¬ 
ary condition at the plate into the Green function -30] 
(Fig. 1ST]) . We assumed a distance between the particle 
tip and the wall of 0.1 pm and checked that the choice of 
this value only has a minor influence on the calculated 
coefficients. In the calculation we restricted the motion 
to the direction along the x-axis in which active oscilla¬ 
tions take place. The resulting coupling coefficients are 
shown in Fig. IS2l In order to reproduce a realistic situ¬ 
ation we simulated two oscillators with a frequency mis¬ 
match H 2 /Hi = 1.2 (based on the measured average fre¬ 
quency difference between particles that are out of range 
of hydrodynamics interactions) and additional noise. The 
resulting correlation coefficients are shown in Fig. 0^. 
They show qualitative agreement with the measured val¬ 
ues under similar conditions, both in their distance- and 
direction-dependence. Frequency mismatch and noise 
both contribute to the loss of synchrony at larger dis¬ 
tances. The relative weight of different coupling coeffi¬ 
cients ctt etc. depends sensitively on the parameter e. 
A small variation can lead to a regime where the two os¬ 
cillators show anti-phase synchronization (Fig. 0)i). As 
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FIG. 5: a) Particles arranged in a row. b) Phase profile: laser 
power 40 mW per particle and distance between adjacent par¬ 
ticles 8-9 pm. c) Simulation of the same system using param¬ 
eters a = 0.8, 0 = 1.5, e = 12 pm. d) Correlation matrix 
of experimentally measured particle positions, e) Equivalent 
correlation matrix from simulation. 


one can see from the phase diagram for typical param¬ 
eters (Fig. [3j dashed line) the parameter range that al¬ 
lows anti-phase synchronization is rather narrow. This 
provides a possible explanation for the rare instances in 
which anti-phase synchronization was observed in the ex¬ 
periment. 

Synchronization in particle row. Following the two- 
particle experiments we investigated if hydrodynamic in¬ 
teractions can also lead to synchronization in a row of 
particles. When biological cilia are densely covering a 
surface their oscillations usually form metachronal waves 
whose direction and wavelength vary from system to sys¬ 
tem. A problem that arises in theoretical models of wave 
formation is that the cilia at the ends of a chain are sub¬ 
ject to different interactions with their neighbors than 
those in the middle. As a consequence some models 
obtain synchronization or metachronal waves when they 
introduce periodic boundary conditions but not in a fi¬ 
nite chain [31] , although the latter is possible, too [32J. 
Likewise, experimental studies on model systems concen¬ 
trated on particles arranged in closed ring structures [l8|. 

In our experiment 7 optical traps were arranged in a 
linear row with a spacing of 8 pm (Fig. [5jt). Particles 
never showed complete synchrony or phase waves. Fig¬ 
ure EJ> shows an example where a high degree of local 
synchrony persisted over many oscillations but not in¬ 
definitely. The correlations are stronger in the middle of 
the chain and weaker at its ends (Fig. [5]i). 

We also extended the theoretical model to a row of 
particles. We used the same parameters as for two par¬ 
ticles and applied the distance-dependent coupling coef¬ 
ficients (Fig. S2) to describe the interaction between a 
particle and all other particles in the chain. A phase plot 
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of a simulated system with the same parameters we used 
for two particles is shown in Fig. [5fc and the correlation 
matrix between all oscillators in Fig. [5^. Like in the ex¬ 
periment partial order but no complete synchronization 
can be seen. The loss of synchrony is mainly due to ge¬ 
ometric effects in the chain, because the particles at the 
end are subject to weaker interactions than those in the 
middle. 

In conclusion, we have shown that elliptical colloidal 
particles oscillating in a laser beam provide a good model 
system for studying hydrodynamic synchronization. In a 
two-particle system we showed that the coupling strength 
depends on the interparticle distance as well as their 
spatial arrangement. While particles usually synchro¬ 
nize in-phase, anti-phase synchronization was observed 
as well. For many particles arranged in a row we observe 


only weak synchronization in agreement with theoretical 
results. A question that remains to be investigated is 
whether 2-dimensional arrays of oscillators would show 
more robust synchronization and what are the necessary 
conditions for the formation of metachronal waves in such 
a system. 
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Particle 2 



FIG. SI: Boundary element representation of two ellipsoidal 
particles (side view) in close proximity to a wall with no-slip 
boundary condition. 
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FIG. S2: Reduced coupling coefficients ctt, crr , ctr , crt calculated with the boundary element method (BEM) as a function 
of interparticle distance d for longitudinal, diagonal and lateral spatial arrangement. 

























